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It is found that the radius of a stable strangelet decreases as the tem- 
perature increases in a quark mass density-dependent model. To over- 
come this difficulty, we extend this model to a quark mass density- and 
temperature- dependent model in which the vacuum energy density at zero 
baryon density limit B depends on temperature. An ansatz which reads 
B = B [l- a{T/T c ) + b(T/T c ) 2 ] is introduced and the regions for the best 
choice of the parameters are studied. 
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It has been expected that a new state of matter, namely, the quark-gluon plasma (QGP) 
can be formed in the relativistic heavy-ion collision. Many signatures of the formation of 
QGP, such as J/^f suppression, strangeness enhancement, thermal dilepton and electro- 
magnetic radiation etc. have been suggested JIJ. A lot of them have been found in recent 
experiments as pointed by CERN collaboration. However whether QGP has been got is 
still an open question, because many signatures suggested can also be explained by other 
treatments B. To search an unambiguous signature is an essential task for the study of 
QGP. Due to this reason, the strangelets, the strange quark matter (SQM) in finite lumps, 
have attracted much attention in recent years. Greiner and his coworkers |3| argued that 
strangelets might be produced in ultrarelativistic heavy-ion collisions and could serve as an 
unambiguous signature for the formation of QGP. 

Since the speculation of Witten |3[ that the strange quark matter (SQM) would be 
more stable than the normal nuclear matter and nuclei, many models including MIT bag 
model 0-0, quark- meson coupling model ||, quark mass-density dependent model (QMDD) 
PHl5|) Friedberg-Lee model |L6|], chiral SU(3) quark model |L7) and etc. jTj| have been 



employed to predict the behavior of SQM in the bulk. It is of interest to investigate 
whether above models which are successful for bulky strange matter can be used to de- 
scribe strangelets. 

This paper evolves from an attempt to study this problem for the QMDD model. This 
model can provide a dynamical description of the confinement mechanism and explain the 
stability of SQM successfully via the suggestion of a density dependent masses for u, d 
and s quarks. It is found that a difficulty emerges if the QMDD model is used to describe 
strangelets at finite temperature. The radius of strangelet decreases as the temperature 
increases. This is of course unreasonable. Since the mass of hadrons is observed to be 
dependent on temperature, we extend the QMDD model to a quark mass density- and 
temperature-dependent (QMDTD) model, and show that after a suitable choice of the ad- 
justed parameters in the function of B(T), the radius of the strangelet increases with the 
rise of temperature. 

According to the QMDD model, the masses of u, d quarks and strange quarks (and the 
corresponding anti-quarks) are given by P,|10 



m g = ^^> (q = u,d,u,d), (1) 
3n B 

B 

m s> s = m s0 + , (2) 
3n B 

where hs is the baryon number density, m S Q is the current mass of the strange quark and 
B is the vacuum energy density inside the bag. 
The thermodynamic potential is 

Q = Y i n i = -Y l T dk^ In (l + e -/*0*(*)-«)) , (3) 

i i 

where i stands for u, d, s (or u,d,s ) and the electron e(e + ), /ij is the corresponding chemical 
potential. Since the masses of quarks depend on density, the thermodynamic potential Q 
is not only a function of temperature, volume and chemical potential, but also of density. 
How to treat the thermodynamics with density-dependent particle masses self-consistently 
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is a serious problem and has made many wrangles in the references JlO ,|TT , 13|, 1 5 1 . Here after 



we follow the treatment of ref. [13|. The number density n i} the total pressure p and the 
total energy density e are given by 
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respectively. The baryon number density ub reads 



(6) 



n B = -^{n u + n d + fi 8 ). 
The conditions of charge neutrality and chemical equilibrium are 

2h u = h d + n s + 3n e , 



(7) 



(8) 
(9) 



Now we are in a position to use the above formulation to investigate the strangelets. 
Instead of one of plane wave, the density of states — — - for a sphere with radius R is needed 



ea- 



rn our calculations. It was given by a multi-reflection theory |T9 



dk 



9i 



2tt 2 



+ f t 



mi 
k 



kS + £ 



(0 



mi 
k 



C + 



(10) 



where V = ^nR 3 ,S = 4ttR 2 ,C = 8nR,gi = 6 for quarks and antiquarks, — 2 for e and 
e + .The second term on the right hand side of Eq. ([H]) corresponds to the surface contribution. 
It is shown 



k J 8n \ n rrii 



(11) 



This term is zero for massless quarks. The third term on the right hand side of Eq.(10) 
comes from curvature of the bag surface. It can not be obtained by the multi-reflection 
theory directly except for two limiting cases — > and rrii — > oo. Madsen proposed [^TJ 
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This simple form agrees to the above mentioned two values in the corresponding limits. 
For the strangelets, the stability condition reads 
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dF 

mi- ' < 13 > 

where F is the total free energy. And it can be given by 

F = E-TS } (14) 

where E = eV is the total energy, 

poo dN- 

S = YlSi = -22J dk— [rii \ia rii + (1 - rii) \n(l - rii)] (15) 

is the entropy, rii is the distributing function of fermions and % stands for u, d, s (or u, d, s) 
quarks and the electron e(e + ). Substituting Eqs.([lO|), (|TlD , fll2"D into Eq.@ and using 
Eq.(O), we can obtain the stable radius R of the strangelet, which is a function of temper- 



ature. 

The curves for F per baryon number A vs. R of QMDD model at zero temperature and 
at T = 50MeV are shown in Fig.l by solid line and dashed line, respectively. The values of 
A, B, m sQ are chosen as 

A = 20, B = 170MeVfm" 3 , m s0 = 150MeV. (16) 

It is clear that the radius for minimum F decreases as temperature increases. The tem- 
perature dependence on radius R is displayed in Fig. 2, in which it can be seen that R(T) 
is a monotonously decreasing function. This is of course unreasonable because the bag is 
expected to expand as the increase of temperature. 

To understand this result physically, let us recall what happens for the masses of nucleons 
and mesons when temperature increases. In fact, the effective masses of nucleons, effective 
masses and screening masses of mesons, are all dependent on temperature |22|-p4"|| . We can 
sum the tadpole diagrams and the exchange diagrams for nucleons, the vacuum polarization 
diagrams for mesons by Thermo Field Dynamics and find that the masses of nucleons and 
mesons all decrease as temperature increases. This result for p— meson is in agreement 



with recent experiments [23] . According to the constituent quark model, the nucleon is 
constructed by three quarks and the meson by two quarks. This means that to satisfy quark 
model we must consider the temperature dependence of the quark mass. But this effect has 
not been taken into account in Eqs.(|IJ) and (Q) if B is a constant. Therefore we come to a 
conclusion that the unphysical result of QMDD model in studying strangelets comes from 
the assumption that B is temperature independent. 

Although it is possible to find the function B(T) from a calculations of the vacuum 
energy, but in this paper, instead of the vacuum energy calculation, we study this problem 
more generally by introducing an ansatz, namely 



B = B n 



TJ \TJ 



(17) 



where a, b are two adjust parameters, and T c = 170MeV is the critical temperature of the 
quark deconfinement phase transition. The reasons for our choice are: first, the calculation 
of B(T) are mo del- dependent and the results given by different model are very different 
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26] p7| ; secondly, almost all previous calculations are based on MIT bag model, but now 
we discuss the QMDD model. In fact we can imagine that the Eq. ([l7|) is the temperature 
expansion of B in the low temperature regions. Since B is zero when T = T c , a condition 



l- a + b = (18) 

is imposed and only one parameter a can be adjusted. 

Introducing Eqs.([17|) and fllBP , we extend the QMDD model to a quark mass density- and 
temperature- model(QMDTD). The results of our model obtained by substituting Eqs.flTTD 
and (|I8"D into Eqs.(|l|) and (|2|) are shown in Figs.(l), (3) and (4). The F/A vs. stable R curve 
for T = 50MeV, a = 0.65 given by QMDTD model is shown in Fig.(l) by dot line. We see 
that the value of stable radius increases from R(r=o) — 2.27fm to -R( T=50 MeV) = 2.31fm, but 
decreases to -R( T=50 MeV) = 2.18fm for QMDD model as shown by dashed-line. 

The value of parameter a can affect the result significantly. The range of possible a is 
determined by physical constraints. For example, there are at least two physical constraints: 
(1), the stable radius R should increase with temperature; (2), the energy per baryon E/A 
increases with temperature also. To show the importance of first constraint, we plot the 
R(T) curves for a = -0.20, 0.20, 0.40, 0.65, 0.80 in Fig. (3), respectively. We see that R(T) 
becomes a monotonously increasing function in the regions < T < 80MeV when a > 0.65. 
On the other hand, the E/A vs. T curves for a = 1.50, 0.90, 0.80, 0.65, 0.20 are shown 
in Fig. (4) to investigate the relevance of the second constraint. We find E/A vs. T curves 
becomes monotonously increasing function in the same temperature region when a < 0.8. 
Therefore the best values for parameter a are in the range 

0.65 < a < 0.8. (19) 

In summary, in order to overcome the difficulty related to the reduction of R with T, 
we suggested a QMDTD model. An ansatz for the temperature dependence of vacuum 
energy B (Eq.flTTD) was introduced. We found that the parameters a must lies in the range 
0.65 < a < 0.8. Our model can be used to study strangelet, in addition it can also be 
employed to address systematically the properties of SQM in bulk. 

This work was supported in part by the NNSF of China under construct No. 199750 10, 
No. 19875009, and the Foundation of Education Ministry of China. 



5 



REFERENCES 



[1] S. A. Bass, M. Gyulassy, H. Stoecker and W. Greiner, J. Phys. G25, Rl (1999) and 
references herein. 

[2] For example, see Bin Zhang, C. M. Ko, Bao-An Li, Ziwei Lin and Ben-Hao Sa, Phy. 

Rev. C62, 054905 (2000). 
[3] C. Greiner, P. Koch and H. stocker, Phys. Rev. Lett. 58, 1825 (1987), Phys. Rev. D44, 

3517 (1991). 
[4] E. Witten, Phys. Rev. D30, 272 (1984). 
[5] E. P. Gilson and R. L. Jaffe, Phys. Rev. Lett. 71, 332 (1993). 
[6] J. Madsen Phys. Rev. D47, 5156 (1993) ; ibid D50, 3328 (1994). 
[7] K. S. Lee and U. Heing, Phys. Rev. D49, 2068 (1993). 

[8] P. Wang, R. K. Su, H. Q. Song and L. L. Zhang, Nucl. Phys. A653, 166 (1999). 

[9] G. N. Fowler, S. Raha and R. M. Weiner, Z. Phys. C9, 271 (1981). 
[10] S. Chakrabarty, Phys. Rev. D43, 627 (1991); ibid D48, 1409 (1993). 
[11] O. G. Benrenuto and G. Lugones, Phy. Rev. D51, 1989 (1995). 
[12] G. X. Peng, P. Z. Ning and H. Q. Chiang, Phy. Rev. C56, 491 (1997). 
[13] G. X. Peng, H. C. Chiang, B. S. Zou, P. Z. Ning and S. J. Luo, Phy. Rev. C62, 025801 

(2000) . 

[14] G. X. Peng, H. C. Chiang, P. Z. Ning, B. S. Zou, Phy. Rev. C59, 3452 (1999). 
[15] P. Wang, Phy. Rev. C62, 015204 (2000). 

[16] Y. J. Zhang, S. Gao, Y. B. He, R. K. Su and X. Q. Li Chinese Phys. Lett. 14, 89 (1997). 
[17] P. Wang, Z. Y. Zhang, Y. W. Yu, R. K. Su and H. Q. Song, Nucl. Phys. A688, 791 

(2001) . 

[18] For a review see C. Greiner and J. Schaffner, Int. J. Mod. Phys. E5, 239 (1996). 

[19] R. Balian and C. Block, Ann. Phys. 60, 401 (1970). 

[20] E. P. Gilson and R. L. Jaffe, Phys. Rev. D71, 332 (1993). 

[21] J. Madsen, Phy. Rev. D50, 3328 (1994). 

[22] Y. J. Zhang, S. Gao and R. K. Su, Phys. Rev. C56, 3336 (1997). 

[23] S. Gao, Y. J. Zhang and R. K. Su, Phys. Rev. C53, 1098 (1996). 

[24] P. Wang, Z. Chong, R. K. Su, P. K. N. Yu, Phys. Rev. C59, 928 (1999). 

[25] G. J. Lolos et al, Phys. Rev. Lett. 80, 241 (1998). 

[26] M. C. Birse, J. J. Rehr and L. Wilets, Phys. Rev. C38, 359 (1988). 

[27] R. D. Pisarski, Phys. Lett. B110, 155 (1982). 

I. FIGURE CAPTIONS 

Fig.l The total free energy per baryon F/A as a function of radius R, for temperature 
T = 0(solid line), T = 50MeV in QMDD model (dashed line) and T = 50MeV, a = 0.65 in 
QMDTD model (dotted line). 

Fig. 2 The radius R as a function of temperature T for QMDD model. 

Fig. 3 The radius R as a function of temperature T for QMDTD model with various 
values for the parameter a as indicated. 

Fig. 4 The energy per baryon E/A as a function of temperature T for QMDTD model 
with various values for the parameter a as indicated. 
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